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We solve the Faddeev-Yakubovsky equations for 3N and 4N bound states based on 
the most modern reaUstic nucleon-nucleon interactions. We include different realistic 3N 
forces. It is shown that all 3N force models can remove the underbinding of the triton and 
a-particle which one obtains with existing NN interactions. The agreement of theoretical 
predictions and the experimental binding energy is quite good and there is little room left 
for the action of four-nucleon forces in the a-particle. The effect of 3N forces on the wave 
function is investigated. 

1. Introduction 

In recent years it has become possible to solve rigorously the 4N Faddeev-Yakubovsky 
equations in momentum space Using parallel computers one can achieve converged 

solutions for realistic nuclear Hamiltonians. There are several realistic interactions avail- 
able which can describe the nuclear two-body problem with a very high accuracy. These 
are the potentials of the Nijmegen group Nijm I, II and 93 0], the CD-Bonn and 
the AV18 interaction f^. All these interactions account for charge independence break- 
ing (CIB). The CD-Bonn and AV18 additionally distinguish pp and nn forces. It is well 
known that nonrelativistic calculations based on these potentials lead to an underbinding 
of 500-800 keV in the 3N system. Two additional dynamical ingredients should cure this 
underbinding: relativistic corrections and three- nucleon forces (3NF). In this contribution 
we would like to address the question whether 3NF's can solve the underbinding problem 
in both the 3N and 4N system and investigate the dependence of binding energies and 
wave functions on the NN and 3N force model. 



2. Binding energies 

In Table |1] we compare our theoretical predictions for various nuclear potentials to 
the experimental binding energies. The results depend on the NN force model, but the 
experimental binding energy of the triton is always underpredicted by 500-800 keV for all 
nuclear potentials. The weakest binding is predicted by the three local potentials AV18, 
Nijm 93 and Nijm II, whereas the strongly non local CD-Bonn gives ~ 300 keV more 
binding. The non local potentials are less repulsive for small distances and their short 
range correlations are softer. This is visible in their smaller kinetic energy predictions P|. 
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Table 1 

3N and 4N binding energies for various NN potentials with expectation values T of the kinetic energy. 





3H 


^He 


Potential 


EB[Me\] 


r[MeV] 


£;B[MeV] 


r[MeV] 


CD-Bonn 


-8.012 


37.42 


-26.26 


77.15 


AV18 


-7.623 


46.73 


-24.28 


97.83 


Nijm I 


-7.736 


40.73 


-24.98 


84.19 


Nijm II 


-7.654 


47.51 


-24.56 


100.31 


Nijm 93 


-7.661 


45.60 


-24.53 


95.34 


exp. 


-8.48 


-28.30 



Table 2 

Cutofl parameters A, 3N and 4N binding energies for various force combinations with expectation values 
T of the kinetic energy. 







3h 


^He 


Potential 


A[m^] 


£;s[MeV] 


T[MeV] 




T[MeV] 


CD-Bonn+TM 


4.784 


-8.480 


39.10 


-29.15 


83.92 


AV18+TM 


5.156 


-8.476 


50.76 


-28.84 


111.84 


AV18+TM' 


4.756 


-8.444 


50.55 


-28.36 


110.14 


AV18+Urb-IX 




-8.478 


51.28 


-28.50 


113.21 


Nijm I+TM 


5.035 


-8.392 


43.35 


-28.60 


93.58 


Nijm II+TM 


4.975 


-8.386 


51.02 


-28.54 


113.09 


exp. 




-8.48 


-28.30 



A linear correlation found in 0] for simple forces holds also in case of the modern realistic 
forces including charge independence and symmetry breaking 0. 

The inclusion of 3NF's cannot result in parameter free predictions for the 3N and 
4N system because all realistic potentials are not accompanied by theoretically founded 
3NF's based on the same fundamental dynamical assumptions. Therefore we use the 
Tucson-Melbourne(TM) 27r 3NF and adjust the cutoff parameter in the NNvr form 
factors individually for each NN interaction insuring that the resulting 3N Hamiltonians 
reproduce the experimental triton binding energy. The Nijmegen forces do not specify a 
neutron-neutron force. In this case we used ^He for the adjustment. The resulting cutoff 
values are given in Table || . In ||10| a modified version of the TM force was suggested 



which at least does not violate chiral symmetry. We call it TM' and adjusted it to the 
AV18 force in the same manner. Additionally, we included the Urbana-IX 3N interaction 

As shown in Table |] all force combinations overpredict the a-particle binding energy. 
However the overbinding is rather small. This indicates that the contribution of 4N forces 
is small in the a-particle. Of course, on this phenomenological level the addition of more 
structures to the 3NF (for example the vr — p exchange) might change this situation. 
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Table 3 

Expectation values of the TM, TM' and Urbana-IX force calculated with wave functions we generated 
using the AV18, AV18+TM,AV18+TM' and AVlS+UrbanalX force combinations. The cutoff values of 
TM and TM' are adjusted to the AVIS potential and given in Table |. 





<TM> 


<TM'> 


<Urb-IX> 


^'(AVIS) 




-1.60 


-3.13 


-2.64 


*(AV18^ 


-TM) 


-8.03 


-4.76 


-4.20 


^'(AV18^ 


-TM') 


-2.99 


-5.17 


-4.98 


*(AV18^ 


-Urb-IX) 


-3.31 


-5.48 


-5.94 



3. Wave function properties 

In this section we address the question whether 3NF's have got any effects on the wave 
functions. Because the effects in the 3N and 4N systems are very similar and are more 
pronounced in the 4N system, we restrict ourselves to the a-particle. 



It is well known that 3NF's are not treatable by perturbation theory ||T2|,|T3|. Therefore 
the inclusion of 3NF's change the wave function. Surprisingly for us we found no significant 
differences in the nuclear momentum distributions and NN correlations comparing wave 
functions generated with and without 3NF. This holds for all force combinations we used. 

To reveal the differences in the wave functions, we calculated expectation values of 
the 3NF operators. In Table |^ we compare the expectation values of the different 3NF's 
obtained with wave functions which we generated using the AV18, AV18+TM, AV18+TM' 
and AVlS+Urbana IX interactions. The results of the wave function generated from the 
AV18 NN interaction alone are generally much smaller in magnitude than the values 
obtained from the "full" wave functions. Clearly the expectation values differ from the 
gain in binding energy one obtains with 3NF's (see Table Q and ||). Therefore the 3NF 
operators are sensitive to the differences in the wave functions. 

In the next step we restrict the wave functions to their S,P and D-wave components 
and calculate the transition matrix elements of the Urbana IX and TM 3NF. The most 
important matrix elements are summarized in Table The contribution of the different 
partial waves is very 3NF dependent. Whereas the repulsive part of the Urbana interaction 
is visible in the diagonal matrix elements, the TM force is strongly attractive in the S- 
S matrix element. The usage of the wave function generated with TM increases the 
attraction in the S-S and S-P transition. Most of the attraction in the Urbana force is 
generated by the S-D transition. In contrast these particular partial waves are extremely 
wave function dependent for the TM force. The results indicate that though the binding 
energies are similar both 3NF operators act very differently on the wave functions of the 
a-particle. 

4. Conclusions 

We presented 3N and 4N binding energies using the most modern nuclear interactions. 
We showed that the underbinding of few-nucleon bound states, which is predicted by 
purely two-nucleon interactions, is solved adding a 3NF. This holds for all combinations 
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Table 4 

Contribution of different parts of the wave function to the Urbana-IX and TM expectation values. In the 
first row the 3N interaction is chosen. In the second row the Hamiltonian, which is used to generate the 
wave function, is fixed and in the first column the part of the wave function on the right and left is given. 
Transition matrix elements are multiplied by 2. The cutoff of TM is adjusted to the AV18 potential and 
given in Table |[ 





Urbana IX 


TM 




*(AV18) ^'(AVISH 


-UrblX) 


^(AVIS) 


^'(AV18- 


FTM) 


S-S 


+3.16 


+2.74 


-2.34 




-4.09 


S-P 


-0.96 


-2.10 


-1.22 




-3.56 


S-D 


-5.44 


-7.46 


+2.08 




-0.14 


D-D 


+0.59 


+0.85 


+0.01 




+0.06 



of 3NF and NN interaction and indicates that 4N forces, at least in the a-particle, are 
small. A look to NN correlations in the wave functions generated from different nuclear 
Hamiltonians reveals that all investigated 3NF's do not change the short range correlations 
in the 3N and 4N system. Because we found at the same time that the expectation values 
of different 3N interactions are strongly dependent on the used wave function, we conclude 
that these wave functions have got very different 3N correlations. The next task of our 
study is to pin down these differences and visualize them. This might help to clarify the 
effects of 3NF's in the 3N and 4N system. 

This work was supported financially by the Deutsche Forschungsgemeinschaft (A.N. 
and H.K.). The numerical calculations were performed on the Cray T3E of the NIC in 
Jiilich. 
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